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Localized  Impact  Damage 
in  Zinc  Sulfide 

by 

H.  P.  Kirchner 
T.  J.  Larchuk 


Ceramic  Finishing  Company 
P.  0.  Box  498 
State  College,  PA  16801 


abstract 


Zinc  sulfide  specimens  were  damaged  by  static  and  impact  loading 
using  glass  spheres.  The  various  types  of  cracks  were  identified  and 
their  sizes  were  measured.  The  maximum  impact  loads  were  estimated  by 
several  methods  including  an  experimental  method  based  on  the  assumption 
that  static  and  impact  loads  are  equal  at  equal  indentation  radii,  a 
semi-empirical  method  previously  developed  on  this  program  and  a  method 
suggested  by  Hutchings.  The  load  (P)  dependence  of  the  crack  length  (c) 


for  static  and  impact  loading  was  best  represented  by 


where  d.  is  the  diameter  of  the  indentation. 
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I.  INTRODUCTION 


Several  types  of  damage  are  induced  when  ceramic  surfaces  are  con¬ 
tacted  by  small  particles.  These  damage  types  include  larger  scale 
damage  such  as  circumferential  ring,  Hertzian  cone,  median,  radial  and 
lateral  cracks,  and  more  localized  damage  such  as  plastic  deformation, 
shear  cracking  and  crushing, d irectly  under  the  contact.  Larger  and 
blunter  particles  favor  elastic  responses  such  as  circumferential  and 
Hertzian  cone  cracks  and  crushing.  Smaller  and  sharper  particles  favor 
elastic-plastic  responses  such  as  median,  radial  and  lateral  cracks, 
and  plastic  deformation  at  the  impact  site. 

Much  of  what  is  known  about  the  mechanisms  of  damage  initiation  and 
propagation  has  been  learned  by  study  of  indentations  formed  under 
static  loading  (hardness)  test  conditions  in  which  an  indenter  is 
pressed  into  a  flat  surface.  Using  small  spherical  tungsten  carbide 
indenters  and  low  loading  rates,  Evans  and  Wilshaw^^  observed  the 
following  sequence  of  crack  formation  in  ZnS: (l)  shallow  radial  surface 
cracks,  (2)  embedded  penny-shaped  cracks  perpendicular  to  the  surface 
(median  vents),  (3)  merging  of  the  median  vents  and  radial  cracks  under 

coplanar  conditions,  and  (4)  subsurface  cracks  approximately  parallel 

(2) 

to  the  surface  (lateral  vents).  Subsequently,  Richard  and  Kirchner 

used  results  of  static  indentations  and  impacts  of  larger,  3  mm  diameter, 

glass  spheres  on  ZnS  for  a  preliminary  evaluation  of  a  theory  describing 

the  elastic-plastic  response  to  impact  damage  and  Shockey,  Dao  and 
(3) 

Curran  investigated  nucleation  and  growth  of  cracks  at  static  inden¬ 
tations  and  impacts  of  small  tungsten  carbide  spheres  on  ZnS. 


Several  mathematical  models  are  available  that  can  be  used  to  predict 


u 


or  scale  the  extent  of  damage  due  to  static  indentation  or  localized  im¬ 
pact  (Table  1)^  Use  of  these  models  depends  on  knowledge  of  the 

contact  load.  Direct  measurement  of  the  contact  load  is  relatively  easy 
for  static  indentations  but,  in  the  impact  case,  it  is  difficult  at 
best. 

In  the  present  investigation,  ZnS  specimens  were  damaged  by  static 
and  impact  loading  using  glass  spheres.  The  Young's  modulus  of  glass 
is  only  about  70/.  as  great  as  that  of  ZnS,  so  the  glass  is  expected  to 
be  subject  to  substantial  elastic  deformation.  On  the  other  hand,  the 
yield  stress  of  ZnS  is  not  very  high  compared  with  other  ceramics,  so  an 
elastic-plastic  response  is  expected.  The  contact  damage  was  character¬ 
ized.  The  impact  loads  were  estimated  using  an  experimental  method 
based  on  the  assumption  that  the  static  and  impact  loads  are  equal  at 

equal  indentation  radii,  and  analytical  methods  using  equations  derived 

(2)  (12) 

by  Richard  and  Kirchner  and  Hutchings  .  The  variations  in  the 
extent  of  damage  were  scaled  using  selected  models  from  Table  I  and  the 
results  were  compared  with  experimental  data. 


^pupijg.  rr 
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TABLE  X  Crack  Size  Relations  for  Various  Contact  Conditions 


TABLE  I  continued 
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II.  PROCEDURES 

Zinc  sulfide  plates  (25  x  25  x  7  mm.)  were  indented  by  glass 

•irk 

spheres  ,  3  mm.  diameter,  under  various  loads  using  an  Instron  testing 
machine.  Before  testing,  the  plates  were  coated  with  a  very  light  layer 
of  soot.  The  contact  boundary  was  marked  by  the  impression  left  in  the 
soot  layer  when  the  sphere  was  removed.  The  contact  radii,  indentation 
radii,  numbers  of  radial  cracks,  and  maximum  radial  and  lateral  crack 
lengths  were  measured  using  a  calibrated  eyepiece  in  an  optical  micro¬ 
scope,  for  a  load  range  of  84  to  1059  N.  Measurements  of  the  maximum 
radial  crack  depth  were  made  by  viewing  through  the  edge  of  the  plates. 
Some  test  sites  were  too  far  from  the  edge  to  be  measured. 

Similar  plates  were  impacted  by  glass  spheres  accelerated  by  a  gas 
gun  using  compressed  ^  gas.  The  impact  and  rebounding  velocities  were 
measured  using  photographs  with  multiple  stroboscopic  images  of  the  glass 
spheres.  The  spheres  were  coated  with  a  thin  layer  of  aluminum  to  im¬ 
prove  their  reflectivity.  The  impact  velocity  range  was  19  to  65  m* **s-^. 
The  impact  damage  was  characterized  as  described  above  for  the  static 
indentations.  In  some  cases  the  profiles  of  the  static  and  impact 
indentations  were  compared  using  a  prof ilometer.  The  results  were 
analyzed  and  are  presented  in  the  following  section. 


*Chemical  vapor  deposited  zinc  sulfide  from  Raytheon  Company, 
Waltham,  MA. 

**No.  3000,  Walter  Stern  Inc.,  Port  Washington,  NY. 
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III.  RESULTS  AND  DISCUSSION 


A.  Comparisons  of  static  and  impact  loading  damage 

The  response  of  the  ZnS  plates  to  static  and  impact  loading  was 
elastic-plastic,  yielding  permanent  indentations  with  the  indentation 
radii  increasing  with  load.  The  contact  radii  were  slightly  greater 
than  the  indentation  radii,  as  expected.  Ring  cracks,  radial  cracks 
and  lateral  cracks  were  observed  for  both  static  and  impact  loading 
(Figure  1).  Indentations  were  first  observed  at  static  loads  of  40  N, 
but  a  load  of  at  least  72  N  was  necessary  to  obtain  a  measurable 
indentation.  At  static  indentations  the  load  thresholds  were  84  N  for 
radial  crack  formation  and  899  N  for  lateral  crack  formation.  In  the 
impact  tests  measurable  indentations,  radial  cracking,  and  evidence  of 
lateral  cracking  were  observed  at  the  lowest  velocity  tested,  19  m*s  ^ , 
and  at  all  higher  velocities. 

Damage  at  static  and  impact  sites  was  compared  at  approximately 
equal  indentation  radii  as  shown  in  Figure  1.  When  the  indentation 
radii  are  relatively  small,  as  in  Figures  1  A  (static)  and  1  B  (impact), 
the  radial  cracks  formed  at  the  impact  sites  are  fewer  in  number  and 
smaller  than  those  formed  at  the  static  sites.  However,  in  the  im¬ 
pact  case,  lateral  cracks  initiate  at  smaller  indentation  radii  as 
indicated  by  the  lighter  area  to  the  left  of  the  indentation.  Also,  the 
ring  cracking  is  much  more  visible  at  the  impact  site.  At  larger  inden¬ 
tation  radii,  as  shown  in  Figure  1  C  (static)  and  1  D  (impact),  both  the 
radial  cracks  and  the  lateral  cracks  are  larger  at  the  impact  site  than 
at  the  comparable  static  site.  The  lateral  cracking  at  the  impact  site 
has  produced  chipping  to  the  right  of  the  indentation, while  the  static 
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site  has  no  chipping.  Also,  the  ring  cracking  within  the  impact  inden¬ 
tations  has  become  much  more  pronounced,  causing  the  contour  of  the  im¬ 
pact  indentation  to  appear  somewhat  irregular,  while  the  static  indenta¬ 
tion  appears  to  be  much  smoother. 

The  irregularities  observed  in  the  impact  indentations  were  investi¬ 
gated  by  profilometer  measurements.  Figure  2  shows  the  diametral  pro¬ 
files  of  static  and  impact  indentations,  compared  at  equal  indentation 
radii.  It  should  be  noted  that  the  vertical  scale  is  ten  times  the 
horizontal  scale.  The  distinguishing  features  are: 

1.  The  impact  indentations  are  much  shallower  than  the  static 
indentations . 

2.  The  centers  of  the  impact  indentations  are  separated  from 
the  outer  rims  by  substantial  circumferential  cracks,  one 
of  which  is  visible  at  the  boundary  of  the  dark  central 
area  in  Figure  1  D. 

3.  The  central  areas  of  the  impact  indentations  have  been 
pushed  down,  apparently  as  a  result  of  shear  failures  at 
the  circumferential  cracks. 

4.  The  outer  portions  of  the  impact  profiles  form  rims  with 
much  lower  slopes  than  those  of  comparable  static  inden¬ 
tations. 

The  features  described  above  indicate  that,  mechanically,  the  forma¬ 
tion  of  impact  indentations  is  much  more  complex  than  the  formation  of 
static  indentations.  These  complexities  are  not  presently  understood. 

It  is  likely  that  some  of  these  differences  occur  because  of  the  strain 

(13 ) 

rate  dependence  of  the  flow  stress  .  Because  the  flow  stresses  are 
higher  during  the  early  stages  of  impact  when  the  strain  rates  are  high, 


11 


StOtiC  lilt 
r-  *  413pm 

Load  *  l,059N 

A  =  35.5  pm 


Impact  site 
Tj  *  413  pm 
V0  *  65.0  m/s 

A  =  23.9  pm 


Static  site 
rj  =  400  pm 
Load  *  936N 
A  =  30.4  pm 

Scales: 


Impact  site 
fj  s'  396  pm 
V0  =  61.6  m/s 
A  =  21.6  pm 

Note  difference 
between  horizontal 
and  vertical 
scales. 


50pm 


i - - 

500  pm 


Figure  2  Diametral  profilometer  traces  of  static  and  impact  indentations 
(ZnS,  3  mm  diameter  glasB  spheres). 
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the  response  is  more  elastic  than  it  is  in  the  static  case.  In  metals, 
in  which  the  yield  stresses  are  a  small  fraction  of  Young's  modulus,  one 
might  expect  the  response  to  be  primarily  plastic  despite  the  increase  in 
flow  stress.  However,  in  ceramics  in  which  the  yield  stresses  are 
usually  a  much  larger  fraction  of  Young's  modulus,  the  response  may  be 
primarily  elastic  during  the  early  part  of  the  impact  when  strain  rates 
are  high.  This  may  account  for  the  reduced  volume  of  the  indentations 
and  the  smaller  indentation  depths  in  the  impact  specimens.  Also,  one 
would  expect  this  greater  resistance  to  deformation  of  the  ZnS  to  be 
accompanied  by  greater  elastic  deformation  of  the  glass  spheres  at  each 
stage  of  the  Indentation. 

One  can  speculate  that  the  boundary  formed  by  the  wide  circumfer¬ 
ential  crack  represents  a  boundary  between  a  central  region  in  which 
forces  tending  to  cause  radial  flow  of  the  ZnS  along  the  interface 
between  the  ZnS  and  glass  are  insufficient  to  overcome  the  frictional 
forces and  a  circumferential  band  in  which  the  frictional  forces  are 
overcome  so  that  movement  of  ZnS  radially  along  this  interface  does 
occur. 

The  existence  of  the  deeper  indentations  at  the  static  sites  implies 

that  the  residual  loads  acting  at  these  sites  after  the  applied  load  is 

removed  are  larger  than  those  at  the  impact  sites.  Such  a  difference  in 

residual  stresses  may  lead  to  observable  differences  in  radial  crack 
(7) 

propagation  . 

The  sketches  in  Figure  3  present  details  of  cracks  observed  by 
examining  cross-sections  of  static  and  impact  loaded  specimens.  The 
pairs  of  specimens  have  comparable  indentation  radii  and  are  arranged  in 
order  of  increasing  indentation  radius.  The  cross-section  in  Figure  3  A 
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Static  site 


Static  site 


Load  *  296N 
fj  *  247  pm 


Load  *  708N 
rj  *  348  pm 


A. 


C. 


Impact  site  Impact  site 


V0  *  19.7  m/s 
fj  *  243pm 


V0  =  40.8  m/s 
r  j  =  351  pm 


B. 


D. 


Figure  3 


Comparisons  of  cracks  formed  at  static  and  Impact  sites  with  comparable 
Indentation  radii  (ZnS,  3  mm  diameter  glass  spheres). 
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represents  damage  at  a  small  static  indentation.  The  damage  consists  of 
radial  cracks  extending  from  the  edges  of  the  indentation  and  a  slightly 
porous  region  under  the  indentation.  The  radial  cracks  are  not  neces¬ 
sarily  coplanar  and  they  do  not  bridge  the  gap  under  the  indentation. 

Evans  and  Wilshaw^^  observed  that  radial  cracks  form  on  loading. 

This  observation  is  confirmed  by  observation  of  radial  cracks  that  extend 
into  the  indentation  (Figure  1  A).  If  the  radial  cracks  formed  on  unload¬ 
ing,  as  a  response  to  the  residual  stresses,  the  radial  cracks  would  be 
expected  to  originate  at  the  boundary  of  the  irreversibly  deformed  zone. 

Formation  of  radial  cracks  on  loading  implies  that  the  out-of -plane 

hoop  stresses  (Figure  4)  are  tensile  even  when  the  indentations  are  very 

shallow.  This  condition  conflicts  with  the  conditions  indicated  for 

Boussinesq  (point  load  at  surface)  and  Hertzian  (blunt  load  at  surface) 

stress  distributions ^ ^  and  the  results  of  calculations  based  on  the 

(9) 

Mitchell  equations  which  indicate  that  the  out-of -plane  hoop  stresses 
are  compressive  near  the  surface.  For  an  extreme  case  in  which  a  rigid 
sphere  is  pressed  to  the  depth  of  its  radius  into  the  surface,  numerical 
analyses  of  the  elastic  stress  field  by  Evans  and  Wilshaw^^  have  shown 
that  the  out-of-plane  hoop  stresses  are  tensile.  If  one  considers  that 
the  irreversibly  deformed  zone  acts  in  some  ways  like  this  embedded 
sphere  to  transfer  the  load  into  the  interior  of  the  specimen,  it  is 
understandable  that  the  out-of-plane  hoop  stresses  can  be  tensile  even 
though  the  indentations  are  shallow.  Friction  may  be  another  contribut¬ 
ing  factor. 

The  radial  cracks  initiate  at  the  surface  and  spread  outward  along 
the  surface  but  they  do  not  spread  downward  into  the  material,  appreci¬ 
ably,  until  they  have  grown  out  past  the  final  indentation  radius. 


Load  axis 


The  surface  of  the  plate  is  represented  by  the  y-z  plane.  Due  to 
symmetry,  the  reference  plane  is  an  arbitrary  plane  perpendicular  to 
the  y-z  plane  and  containing  the  loading  axis  (x-axis).  Then,  at  point 
o,  ar  is  the  radial  stress,  <Jq  is  the  in-plane  hoop  stress  and  is 
the  out-of-plane  hoop  stress. 


Figure  4  Coordinate  axes  and  identification  of  stresses. 
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The  comparable  impact  site  (Figure  3  B)  has  a  somewhat  similar 
radial  crack  on  one  side  of  the  indentation  and  a  small  lateral  crack 
on  the  other.  The  fact  that,  at  small  impact  indentations,  the  radial 
cracks  are  smaller  and  less  numerous  is  consistent  with  the  earlier 
observation  that  the  impact  indentations  are  shallower,  which  would 
reduce  the  out-of-plane  tensile  stresses  during  loading  and  the  residual 
stresses  after  unloading.  The  presence  of  the  small  lateral  crack  at 
the  Impact  site  and  not  at  the  comparable  static  site  which  is  expected 
to  have  higher  residual  stresses,  tends  to  refute  the  argument  that  the 
lateral  cracks  form  in  response  to  residual  stresses.  Instead,  it  is 
likely  that  this  crack  formed  because  the  impact  load  is  greater  than 
the  static  load  for  a  given  indentation  size  so  that  the  lateral  crack 
formed  in  response  to  the  greater  in-plane  hoop  stress.  The  clear 
zone  directly  under  the  indentation  is  shallower  in  the  impact  case  than 
it  is  in  the  static  case.  If  this  clear  zone  is  considered  to  be  a 
plastically  deformed  zone,  perhaps  formed  in  response  to  shear  stresses, 
tne  decreased  depth  is  consistent  with  higher  flow  stresses  at  higher 
strain  rates.  This  observation  is  also  consistent  with  Evans'  obser¬ 
vation  that  lateral  cracking  is  more  extensive  in  harder  materials 
based  on  the  increase  in  dynamic  hardness  with  increasing  strain  rate. 

At  the  next  larger  indentations  (Figures  3  C,  D),  the  radial  cracks 
at  the  static  indentation  are  much  larger  than  those  in  Figure  3  A,  and 
they  have  extended  under  the  indentation  so  that  they  almost  close  the 
gap.  Otherwise,  the  damage  is  similar  to  that  in  Figure  1  A.  The  im¬ 
pact  damage  shows  larger  radial  and  lateral  cracks,  as  expected.  In 
addition,  under  the  clear  zone,  there  are  cracks  extending  deeper  into 
the  ZnS  perpendicular  to  the  surface  that  have  been  labelled  as  median 
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cracks.  These  cracks  are  not  coplanar  and  the  vertical  line  extending 
along  the  impact  axis  represents  the  intersection  of  the  cracks.  Median 
cracks  are  frequently  observed  at  static  indentations  under  sharp  inden- 
tors.  As  the  load  is  removed,  radial  cracks  may  extend  in  response  to 
residual  stresses  and  merge  with  existing  median  cracks.  Apparently, 
this  has  occurred  in  the  present  case. 

At  the  largest  indentations  (Figure  3  E,  F),  the  radial  cracks  at 
the  static  indentation  are  even  larger  and  they  have  extended  under  the 
indentation  to  completely  close  the  gap  to  form  an  approximately  semi¬ 
circular  crack.  A  large  lateral  crack  was  observed  close  to  the  sur¬ 
face.  The  impact  damage  shows  larger  non-co planar  radial  cracks  and 
larger  lateral  cracks  compared  with  those  in  Figure  3  D.  Again,  the 
median  cracks  have  merged  with  the  radial  cracks.  Lateral  cracks  are 
observed  at  several  levels,  all  farther  from  the  surface  than  the  lateral 
crack  at  the  comparable  static  site. 

The  numbers  and  lengths  of  cracks  in  static  and  impact  lo  ided 
specimens  are  compared  in  Figures  5  and  6.  Radial  cracks  are  observed 
at  smaller  indentation  radii  at  static  indentations,  and  larger  numbers 
of  radial  cracks  are  formed  at  a  given  indentation  radius  (Figure  5). 

The  number  of  radial  cracks  depends  on  the  total  tensile  strain  at  the 
periphery  of  the  indentation.  Therefore,  observation  of  larger  numbers 
of  cracks  at  static  indentations  indicates  more  total  tensile  strain  in 
these  cases.  Apparently,  the  high  flow  stress  at  high  strain  rates 
reduces  the  total  strain  in  the  impact  cases. 

Radial  cracks  at  static  indentations  are  longer  than  those  at 
impact  sites  for  indentation  radii  up  to  about  400  ^m  (Figure  6).  At 
greater  radii  the  radial  cracks  at  impact  sites  are  longer.  If  the 
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Number  of  radial  cracks  vs.  indentation  radius  for  static  and 
impact  loading  (ZnS,  3  mm  diameter  glass  spheres). 
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vs.  indentation  radius  (ZnS,  3  mm  diameter  glass  spheres). 


21 


arguments  presented  previously  are  correct,  the  lengths  of  the  radial 
cracks  depend  on  the  loads  and  the  depths  of  the  irreversibly  deformed 
zones  which,  in  turn,  influence  the  magnitude  and  size  of  the  out-of- 
plane  hoop  stresses  and  the  residual  stresses.  Apparently,  at  small 
indentation  radii,  the  deeper  indentations  and  greater  plastic  defor¬ 
mation  observed  in  the  static  sites  (leading  to  greater  residual 
stresses)  result  in  greater  out-of-plane  hoop  stresses  than  those 
present  at  comparable  impact  indentations.  However,  at  large  indenta¬ 
tion  radii,  the  larger  loads  necessary  to  overcome  the  high  flow  stresses 
at  high  strain  rates  become  dominant  and  the  cracks  are  larger  than 
those  at  static  indentations. 

B.  Impact  loads 

To  estimate  the  crack  lengths  using  the  equations  in  Table  I,  it 

is  necessary  to  estimate  the  maximum  impact  loads.  Five  methods  of 

estimating  these  loads  are  described  in  Table  II. 

To  determine  the  impact  loads  by  the  first  method,  the  load(P) 

vs.  indentation  radius  (r.)  data  from  the  static  indentation  tests 

1 

were  plotted  (Figure  7).  The  data  were  fitted  by  a  power  law  equation 
using  the  least  squares  method  yielding 

r.  =  2.1O10-5  P0,43  (1) 

i 

The  fact  that  the  exponent  of  P  is  less  than  0.5  indicates  that  P  is 
increasing  faster  than  the  area  of  the  indentation.  In  other  words, 
strain  hardening  is  occurring.  Similarly,  the  impact  velocity  (Vq)  vs. 
indentation  radius  data  from  the  impact  tests  are  presented  in  Figure  8. 


TABLE  II _ Contact  Load  Equations 
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Static  load  vs.  indentation  radius  (ZnS,  3  mm  diameter 
glass  spheres). 


Impact  velocity  vs.  indentation  radius  (ZnS,  3  mm  diameter 
glass  spheres). 
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The  data  were  fitted  yielding 

r  =  6.38-10"5  V0*45  (2) 

i 

Equations  (1)  and  (2)  were  combined  to  obtain  the  impact  velocity 
dependence  of  the  load  which  was 

P  =  13.25  V  1,05  *  (3) 

o 

assuming  that  the  static  and  impact  loads  are  equal  at  equal  indentation 
radii  (Figure  9).  This  calculated  load  increases  slightly  faster  than 
the  velocity.  In  the  elastic  case  (the  Hertz  theory),  the  theoretical 
value  of  the  exponent  is  1.2.  For  the  perfectly  plastic  case,  Hutchings 
found  that  the  maximum  load  was  proportional  to  the  impact  velocity.  It 
is  reasonable  that  the  above  results  fall  between  these  ideal  cases. 

It  should  be  pointed  out  that  this  method  fails  to  account  for 
inertial  effects  and  the  increase  in  flow  stress  with  strain  rate.  Both 
of  these  effects  would  be  expected  to  increase  the  maximum  impact  load 
over  that  calculated,  using  this  method.  However,  the  method  does 
account  for  strain  hardening  to  the  extent  that  it  has  increased  the 
static  indentation  loads. 


Richard  and  Kirchner  method 

.  (2) 

Richard  and  Kirchner  described  a  semi-empirical  method  of 
determining  the  load  in  impact  experiments.  This  method  is  based  on 
the  empirical  relations 


n 


(4) 


(12) 


and 


rj'  = 


(5) 


Impact 
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Figure  9  Impact  load  ve.  impact  velocity  (The  data  points  were  determined 
by  substituting  impact  r*  values  in  equation  (l)  to  calculate 
impact  load  which  was  plotted  vs.  experimental  impact  velocity; 
ZnS,  3  mm  diameter  glass  spheres.). 


26 


in  which  *j  and  *?'  are  the  pressures  averaged  over  the  projection  of  the 
combined  contact  and  indentation  areas  and  the  indentation  areas  respec¬ 
tively,  r  is  the  contact  radius,  r.  is  the  indentation  radius,  r  is  the 
sphere  radius  and  « ,  k * ,  $  and  are  empirical  constants  to  be  deter¬ 

mined  experimentally.  Figure  10  shows  that  the  logarithmic  indentation 
relationship  is  approximately  linear  and  that  representation  of  the 
data  by  the  above  equations  is  reasonable.  Also,  these  results  can  be 
compared  with  those  of  Evans  and  Wilshaw^^  for  smaller  spheres.  The 
contact  parameters  were  determined  from  the  slopes  and  intercepts  of 

plots  of  this  type  with  the  results  shown  in  Table  III  and  compared  with 

(2) 

similar  results  of  Richard  and  Kirchner  .  The  rather  large  difference 
in  the  value  of  * '  between  the  two  investigations  apparently  arose  as  a 
result  of  differences  in  criteria  for  locating  the  indentation  boundaries 
adopted  by  the  person  making  the  measurement  in  each  case. 


TABLE  III  Contact  Parameters 


%  ' 

GNm'2 

% 
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GNm 

Present 

investigation 

0.34 

3.0 

0.71 

4.9 

Richard 

and  Kirchner 

0.27 

5.5 

0.78 

5.5 

and  used  to  calculate  the  impact  load.  Unlike  equation  (3),  this 
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Figure  10  Static  indentation  pressure  vs.  the  normalized  indentation 
radius,  r^/r0  (ZnS  statically  loaded  by  3  mm  diameter 
glass  spheres). 
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method  accounts  for  the  increases  in  pressure  caused  by  inertial  effects. 
However,  it  does  not  account  for  the  increase  in  load  caused  by  other 
material  related  factors  that  may  contribute  to  the  strain  rate  dependence 
of  the  flow  stress  or  the  decrease  in  load  caused  by  the  stress  wave 
energy  loss. 

(12) 

The  third  method  that  was  used  is  Hutchings  method  which  is 
described  in  Appendix  A  and  compared  with  the  first  two  methods.  The 
remaining  methods  listed  in  Table  II  were  not  investigated. 

The  loads  calculated  using  the  first  two  methods  are  plotted  in 
Figure  11,  showing  that  there  are  substantial  differences  in  the  calcu¬ 
lated  load  depending  on  the  particular  method  used  for  the  calculation. 

It  seems  likely  that  the  loads  obtained  using  equation  (3)  are  under¬ 
estimated  because  the  inertial  effects  and  the  strain  rate  dependence  of 
the  flow  stress  are  not  accounted  for.  The  loads  obtained  using  equation 
(6)  may  be  overestimated  or  underestimated.  The  fact  that  the  contact 
radii  are  overestimated,  (Figure  12)  using  the  theory  of  Richard  and 
Kirchner  (the  contact  radii  were  calculated  using  equation  (6)  of 
reference  (2)),  indicates  that  equation  (6)  might  tend  to  overestimate 
the  loads,  but  the  neglect  of  the  strain  rate  dependence  of  the  flow 
stress  by  equation  (6)  would  tend  to  underestimate  them.  In  fact, 
the  lower  contact  radii  observed  experimentally  may  occur  partly  because 
of  this  strain  rate  dependence  of  the  flow  stress  and  partly  because  of 
the  stress  wave  energy  loss. 

In  evaluating  the  strain  rate  dependence  of  the  flow  stress,  it  is 
well  to  remember  that  the  strain  rate  varies  with  time  during  the 
loading  cycle.  Initially,  the  strain  rate  is  high  because  the  velocity 
is  high  but,  as  the  sphere  slows  down,  the  strain  rate  will  decrease  so 
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Figure  11  Impact  load  calculated  by  two  methode  ve.  impact  velocity  CZnS, 
3  mm  diameter  glaas  spheres). 
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that  near  the  end  of  the  loading  cycle,  where  the  maximum  load  is 
observed,  the  flow  stress  may  approach  the  9tatic,  work  hardened  value. 


C.  Load  dependence  of  crack  length 

Examination  of  Table  I  shows  that  several  expressions  are  available 

relating  load  and  crack  length  for  various  assumed  conditions.  Perhaps 

(4) 

the  simplest  approach  was  illustrated  by  Lawn  and  Fuller  who  used  an 
energy  balance  method  to  show  that,  for  blunt  indenters 


P  «  Kic  C 


3/2 


(7) 


This  relation  is  expected  to  apply  only  in  the  far  field;  that  is,  when 
the  crack  length  is  much  greater  than  the  contact  radius.  Also,  even 
though  it  may  describe  the  situation  when  the  load  is  applied,  it  fails 
to  describe  crack  extension  during  unloading  when  residual  stresses, 
originating  at  the  irreversibly  deformed  zone,  substantially  alter  the 
stress  distributions. 

During  ceramic  machining  experiments  in  which  the  contact  dimension 
was  varied  by  varying  the  sharpness  or  bluntness  of  the  diamond  points, 
it  was  observed  that  the  penetration  of  the  damage  depended  strongly  on 


the  contact  dimension  even  when  the  crack  length  was  several  times 

icula 

(10) 


greater  than  the  contact  dimension^^.  In  this  particular  case,  the 


assumption  of  line  contact  loading,  as  represented  by 

p#c  KicL  c% 


(8) 


in  which  L  is  the  contact  length,  yielded  a  satisfactory  comparison  of 

theory  and  experiment.  This  was  the  case  despite  the  fact  that  the  line 

contact  (L)  was  much  smaller  than  the  radial  crack  length  (c).  Subse- 

(2 ) 

quently,  Conway  and  Kirchner  analysed  propagation  of  median  cracks 
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beneath  an  irreversibly  deformed  (disturbed)  zone  and  found  that 

P  «  K  z  c 2  (9) 

IC  o 

where  zq  is  the  depth  of  the  disturbed  zone.  Comparisons  show  that  in 

many  cases  L  and  z^  are  approximately  equal. 

Predictions  based  on  (8)  and  (9)  can  yield  very  different  results 

from  those  based  on  (7).  Therefore,  it  was  of  interest  to  compare  these 

3/2  J  - 

relations  with  experimental  data.  Plots  of  P  vs.  c  and  P/d^  vs.  c 
were  compared,  where  d^  is  the  indentation  width  which  is  assumed  to 
replace  L  or  z^  in  the  above  equations. 

Damage  at  static  indentations 

A  satisfactory  model  relating  load  and  crack  length  should  provide 
for  the  following: 

1.  A  load  threshold  below  which  cracks  do  not  form. 

2.  A  transition  region  in  which  small  variations  in  load 
cause  large  variations  in  crack  length  as  the  crack 
"pops  in"  and  propagates  in  the  near  field. 

3.  A  wide  range  of  loads  in  which  crack  propagation  is 
governed  by  a  simple  relation  based  on  the  far  field 
stress  distribution. 

The  results  should  be  evaluated  based  on  these  criteria. 

3/2  \ 

The  data  are  plotted  as  P  vs.  c  and  P/d^  vs.  c  where  c  is  the 

maximum  radial  crack  length  measured  along  the  surface  from  the  edge  of 

the  indentation  in  Figures  13  and  14.  The  P  vs.  c  curve  shows  a  load 

threshold.  The  data  points  have  a  much  larger  slope  at  low  loads  than 

at  high  loads,  indicating  that  the  exponent  of  c  is  too  high.  The  near 

field  region  is  indicated  approximately  by  c  <  d  where  d  is  the  contact 

c  c 

diameter.  The  variation  of  crack  length  with  load  in  this  region  is  not 
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much  different  from  the  rest  of  the  curve. 

The  P/d.  vs.  c2  curve  is  much  more  consistent  with  the  listed  criteria. 
1 

It  shows  a  load  threshold,  a  transition  region  where  c  <  d^,  and  an 
essentially  linear  dependence  in  the  far  field.  The  coefficient  of  deter¬ 
mination  for  the  data  in  the  far  field  is  0.98,  indicating  an  excellent 
fit.  These  results  provide  additional  experimental  confirmation  of  the 

earlier  results  obtained  in  ceramic  machining  tests  and  assuming  that 

(9) 

d^  o:  z^,  they  support  the  analysis  of  Conway  and  Kirchner 

i. 

Similarly,  a  P/d^  vs.  c2  curve  is  given  in  Figure  15,  where  c  is 
the  radial  crack  depth  measured  under  the  indentations.  The  results 
show  that  the  relationship  between  these  variables  is  approximately  linear. 

Curves  plotted  using  maximum  crack  lengths  measured  from  the  centers 
of  the  indentations  were  less  consistent  with  the  listed  criteria. 

Damage  at  impact  indentations 

The  variation  of  crack  length  with  load  at  impact  indentations  was 

evaluated  similarly,  using  loads  calculated  by  the  methods  described  in 

3/2  \ 

Section  III  B.  The  results  for  P  *  c  and  P/d.  oc  c  where  c  is  the 

i 

radial  crack  length  measured  from  the  edge  of  the  indentation  are  given 

3  /  2 

in  Figures  16  and  17.  The  P  vs.  c  curves  depart  substantially  from 

straight  lines  showing  that  the  exponent  3/2  is  much  too  large.  The  r/d^ 

% 

vs.  c  curves  are  approximately  straight  lines  that  do  not  pass  through 
the  origin.  Based  on  the  results  for  static  indentations  (Figure  14), 
one  would  expect  the  curve  to  pass  through  the  origin.  The  fact  that  it 
does  not  pass  through  the  origin  may  indicate  that  the  loads  in  the  low 
part  of  the  load  range  are  overestimated.  Also,  at  the  impact  sites  there 
may  be  less  crack  extension  on  unloading  because  of  the  smaller  residual 
stresses  at  these  sites  as  a  result  of  the  smaller  irreversibly  deformed 
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Figure  17  Impact  load  divided  by  indentation  diameter  vs.  radial  crack 

length  measured  from  the  edge  of  the  indentation  and  raised  to 
the  1/2  power,  for  two  methods  of  determining  the  load  (ZnS, 

3  mn  diameter  glass  spheres). 
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No  transition  region  was  observed.  However,  one  might  not  observe 
the  transition  region  because  the  minimum  loads  in  the  impact  tests  are 
greater  than  those  for  which  the  transition  was  observed  in  the  static 
tests. 

The  radial  crack  depths  were  measured  in  several  cases.  The 

measurements  did  not  include  the  depths  of  median  cracks,  if  present. 

In  some  cases,  especially  for  smaller  cracks,  the  measurements  were  not 

made  because  the  cracks  could  not  be  observed  clearly  through  the  edges 

of  the  plates.  A  P/d.  vs.  c2  curve,  plotted  using  loads  calculated  by 

(2) 

the  Richard  and  Kirchner  method,  is  presented  in  Figure  IB  showing  a 

linear  variation  and  passing  approximately  through  the  origin. 

Theoretical  crack  lengths  calculated  using  the  two  methods  of 

estimating  the  loads  are  compared  with  experimental  data  for  various 

impact  velocities  in  Figure  19.  The  crack  lengths  used  to  plot  the 

theoretical  curves  were  calculated  by  calculating  the  loads  and  then 

substituting  these  loads  in  the  statically  determined  relation, 

7  h 

P/d^  =  2.74  x  10  c  ,  which  was  obtained  by  a  least  squares  fit  of  the 

data  in  Figure  14.  The  results  show  that,  at  low  impact  velocities,  crack 

lengths  calculated  using  P  =  13.25  (equation  3)  yield  the  best 

estimates  but  at  higher  velocities  the  data  approach  the  predictions  by 

(2) 

the  Richard  and  Kirchner  method.  It  seems  reasonable  that  equation 
(3)  would  yield  its  best  results  at  low  velocities  because  momentum  and 
energy  effects  which  are  neglected  in  this  case  should  be  less  important 
in  this  velocity  range.  It  also  seems  reasonable  that  the  Richard  and 
Kirchner  method  might  yield  better  results  than  equation  (3)  at  higher 
impact  velocities  because  momentum  and  energy  effects  are  accounted  for 
in  this  derivation.  However,  it  is  not  clear  why  the  Richard  and  Kirchner 
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Figure  18  Impact  load  divided  by  indentation  diameter  vs.  radial 
crack  depth  to  the  1/2  power  (impact  load  predicted 
using  equation  (6);  ZnS,  3  mm  diameter  glass  spheres). 
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,  **** 


Radial  Crack  Length 


Figure  19  Comparison  of  two  theoretical  curves  of  radial  crack  length  vs. 

impact  velocity  with  experimental  data  using  c  measured  from 
the  edge  of  the  indentation.  (ZnS,  3  mm  diameter  glass  spheres). 
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method  overestimates  the  crack  lengths  by  such  a  wide  margin  at  low 
velocities . 

D.  Energy  losses 

The  various  energy  loss  mechanisms  during  elastic-plastic  impact 
on  materials  include  the  indentation  energy,  stress  wave  energy  and  the 
fracture  energy  absorbed  by  formation  of  new  surfaces.  Previous 
research  has  shown  that  the  stress  wave  and  fracture  energy  losses  are  not 
particularly  important  factors  .  Hunter predicted  that  for 

elastic  conditions  the  kinetic  energy  loss  should  increase  as  the  2.6 
power  of  the  impact  velocity  as  a  result  of  stress  wave  energy  losses. 

(19) 

Kirchner  and  Gruver  found  that  for  impacts  of  glass  spheres  on  glass 

plates  at  temperatures  below  the  softening  point,  the  measured  kinetic 
energy  loss  variation  was  very  close  to  Hunter's  prediction.  At  temper¬ 
atures  above  the  softening  point  where  indentation  energy  became  important, 
the  energy  losses  increased  but  the  velocity  dependence  remained  very 
much  the  same. 

The  kinetic  energy  loss  vs.  impact  velocity  data  for  impacts  of 
glass  spheres  on  zinc  sulfide  plates  is  given  in  Figure  20.  The  slope 
is  2.94  at  the  higher  velocities  and  the  indications  are  that  the  slope 
is  even  higher  at  low  velocities.  Therefore,  the  kinetic  energy  losses 
increase  more  rapidly  with  velocity  in  ZnS  than  they  do  in  glass.  This 
more  rapid  increase  may  occur  because  of  broadening  of  the  load  vs. 
penetration  curve  with  increasing  impact  velocity  that  can  occur  as  a 
result  of  the  strain  rate  dependence  of  the  flow  stress. 


Energy  Loss 
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IV.  SUMMARY  AND  CONCLUSIONS 

ZnS  plates  were  statically  indented  under  various  load9  and  impacted 
at  various  velocities  using  3  mm  diameter  glass  spheres.  The  contact 
damage  was  characterized  and  the  damage  at  the  static  and  impact  sites 
was  compared  with  the  results  described  in  the  previous  section. 

Damage  at  the  static  and  impact  indentations  included  the  irreversibly 
deformed  zone  directly  under  the  indentation  and  circumferential  ring, 
radial,  and  lateral  cracks.  The  scale  of  the  damage  at  the  static  inden¬ 
tations  increased  with  increasing  loads  on  the  glass  spheres.  An  over¬ 
all  comparison  of  the  static  and  impact  damage  at  approximately  equal 
indentation  radii,  as  shown  in  Figure  3,  shows  that  the  scale  of  the 
impact  damage  is  approximately  the  same  as  that  of  the  static  indentation 
damage.  If  one  assumes  that  the  sizes  of  the  cracks  depend  primarily  on 
the  maximum  loads  applied  at  the  contact,  this  observation  appears  to 
justify  the  original  assumption  of  equal  loads  at  equal  indentation 
radii  for  static  and  impact  indentations.  However,  the  detailed  obser¬ 
vations  show  that  there  are  substantial  differences  in  the  indentation 
radius  thresholds  for  crack  formation  and  the  numbers  and  sizes  of  the 
cracks  at  static  and  impact  indentations  so  that  other  factors  are  impor¬ 
tant. 

At  equal  indentation  radii,  the  impact  indentations  were  shallower 
than  the  static  indentations  showing  that  less  irreversible  deformation 
had  occurred  at  the  impact  sites.  Because  less  irreversible  deformation 
had  occurred,  the  threshold  indentation  radius  for  formation  of  radial 
cracks  was  shifted  to  higher  indentation  radii  and  fewer  radial  cracks 
were  formed.  Also,  at  smaller  indentation  radii  the  radial  cracks  at 
static  indentations  were  longer  than  those  at  impact  sites.  At  the  largest 
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indentation  radii  the  reverse  was  true. 

When  compared  on  the  same  basis,  lateral  cracks  were  initiated  first 
at  impact  sites  and  were  only  observed  at  the  static  sites  with  the  three 
largest  indentation  radii.  At  these  three  static  sites,  the  lateral 
cracks  were  shorter  than  those  at  the  comparable  impact  sites.  The 
lateral  cracks  are  considered  to  form  mainly  in  response  to  the  in-plane 
hoop  stresses.  Therefore,  the  longer  lateral  cracks  at  the  impact  sites 
indicate  that  the  in-plane  hoop  stresses  are  greater  in  the  impact  case. 

Several  mathematical  models  are  available  that  can  be  used  to  pre¬ 
dict  the  lengths  of  cracks  depending  on  the  load  applied  during  contact. 

T/2 

In  the  case  of  static  loading  two  of  these  models,  P  *  (T  and 
P/d^  «•  c  ,  were  compared  with  experiment  by  plotting  the  data  and  deter¬ 
mining  the  intercepts  and  linearity  of  the  resulting  plots.  P/d^w-  c2 

3/2 

yielded  the  best  results.  Linear  extrapolations  using  P  *  c  greatly 
underestimated  the  crack  lengths  at  high  loads.  The  good  results  using 
P/d.  *  c  indicate  that  the  characteristics  of  the  contact  have  a  sub- 

i 

stantial  influence  on  the  extent  of  crack  propagation.  The  mechanism 
responsible  for  this  influence  remains  uncertain.  The  dimensions  of 
the  contact,  depth  of  the  disturbed  zone,  wedging,  and  residual  stresses 
may  be  important  factors  influencing  the  extent  of  damage. 

In  order  to  make  a  similar  evaluation  in  the  impact  case,  it  was 
necessary  to  estimate  the  loads.  Three  methods  were  used,  each  yield¬ 
ing  substantially  different  results.  However,  in  this  case  also,  plot- 
ting  the  data  using  P/d^  vs.  c  yielded  the  best  results. 

To  predict  strength  degradation,  erosion  rates  and  loss  of  trans¬ 
mission  due  to  reflections  at  lateral  cracks,  it  is  necessary  to  predict 
the  extent  of  crack  propagation  for  various  impact  conditions.  The  pre¬ 
sent  results  show  that  for  impacts  of  3  mm  glass  spheres  on  zinc  sulfide 
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at  velocities  up  to  65  m*s  \  the  characteristics  of  the  contact  zone 
have  an  important  influence  in  the  extent  of  crack  propagation,  even  in 


the  far  field 


LIST  OF  SYMBOLS 


P 

c 


k(0 

H 

X 


0 

a 

i(v) 


-  Maximum  load  perpendicular  to  surface 

-  Crack  length  V 

-  Critical  stress  intensity  factor 

-  Dimensionless  function  of  Poisson’s  ratio  v 

-  Hardness  (pressure) 

-  Proportionality  factor  accounting  for  unknown  geometrical  and 
frictional  factors 

-  X  f or  residual  stress  field  after  removal  of  load 

-  Geometrical  factor  determined  by  inelastically  deformed  zone 

-  Geometrical  factor  determined  by  indenter 

2  2  A 

-  >K v )  (l-i/  )  (1-20  /2tt  where  v  is  Poisson's  ratio 
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-  Contact  radius 

-  Radius  of  impacting  sphere 


-  Young' 8  modulus 

-  Q  =  f  (  */,E) 

-  Impact  velocity 


-  n  = 


*  IT  r 


M  ($+4)  ' 
2  k  ir  r  3 
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where  M  is  mass  of  impacting  sphere 
and  k  and  %  are  contact  parameters 


-  Depth  of  "disturbed"  zone 

-  Contact  length 

-  X  for  line  loading 

-  Residual  load  caused  by  irreversible  deformation  during  loading 

-  Indentation  radius 

-  Indentation  diameter 
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Symbols  (continued) 

-  Contact  diameter 

-  Depth  of  the  indentation 

-  Plastic  energy  (energy  of  plastic  deformation) 


-  Young's  modulus  of  projectile 

-  Young's  modulus  of  plate 

-  Poisson's  ratio  of  the  projectile 

-  Poisson's  ratio  of  the  plate 

-  Contact  pressure 

-  Indentation  pressure 

-  Impact  pressure 

-  Area 

-  Resistance  to  penetration 

-  Exponent  characteristic  of  the  variation  of  the  pressure  of 
contact  with  increasing  load 

-  Resistance  to  indentation 

-  Exponent  related  to  work  hardening 

-  Density  of  projectile 

-  Density  of  plate 
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APPENDIX  A 

COMPARISON  OF  THE  TffllEE  IMPACT  LOAD  EQUATIONS 

A.  Hutchings'  impact  load  equation 

Hutchings^)  derived  his  impact  load  equation  by  assuming  that 
during  contact  the  sphere  remains  rigid  and  the  material  responds  per¬ 
fectly  plastically  (no  elastic  response  and  no  work  hardening).  These 
assumptions  enable  Hutchings  to  use  a  constant  indentation  pressure, 
(7^),  in  the  equation  of  motion  of  the  sphere.  The  solution  of  this 
equation  of  motion  yields  a  relation  for  the  load  on  the  sphere  as  a 
function  of  contact  time.  We  are  interested  in  the  maximum  load,  which 
was  found  to  be 

P  =  4  nr  ^  P..  V  (37j/2p1)2.  (l) 

3  o  1  o  d  1 

By  substituting  for  the  density  of  the  sphere,  (p^),  its  mass  divided 

3 

by  its  volume  (3M/47rr  ),  equation  (l)  becomes 

o 

P  =  (2rr  M  r)  I-5  V  .  (2) 

o  d  o 

B.  Comparison  of  Hutchings'  Equation  with  that  of  Richard  and  Kirchner 

(2) 

Richard  and  Kirchner  also  assumed  that  the  sphere  remained  rigid 
during  contact,  but  in  order  to  account  for  the  elastic-plastic  response 
of  the  material  and  the  effects  of  work  hardening,  they  assumed  that  the 
pressure  was  related  to  the  contact  radius  by 

7  =  «  (r  It  .  (3 ) 

c  o 

Then,  by  equating  the  kinetic  energy  of  the  sphere  to  the  work  done  by 


11  <pu 
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the  sphere  in  reaching  maximum  penetration,  the  impact  load  equation 


was  found  to  be 


2  rM(4+4>" 

;  r  r  - =- 

°  2**r  3 

L  O  J 


(^v  2(«). 


If  the  assumptions  made  by  Hutchings  are  applied  to  equation  (3),  £ 
becomes  zero,  and  K  then  equals  a  constant  indentation  pressure,  7}^. 
Making  these  substitutions  for  $  and  «  ,  equation  (4)  becomes 


4M 

2*d*r  3 
L  ‘a  o  J 


which  reduces  to 


P  =  (2nr  Mi?,)2  V  . 

o  d  o 


Tims,  the  impact  load  equation  of  Richard  and  Kirchner  reduces  to 
Hutchings'  equation  for  the  ideal  case  of  a  rigid  sphere  and  a  perfectly 
plastic  material  response. 


C.  Comparison  of  the  three  impact  load  relations 

To  see  how  the  inclusion  of  elastic-plastic  response  and  work 
hardening  affects  the  impact  load,  equations  (2)  and  (4)  were  evaluated 

(3) 

and  compared.  To  evaluate  equatior  (2),  the  hardness  (H  =  1.90  GPa  ) 
was  substituted  for  Thus  equation  (2)  reduces  to 

P  =  25.2  V  .  (7) 

o 

When  the  experimental  values  for  $  and  k  are  substituted  into  equation 
(4),  we  get 


P  =  16.1  V 
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Equations  (7)  and  (8)  are  plotted  in  Figure  1.  For  velocities  up  to 
about  25  m*s  \  equation  (7)  is  a  good  linear  approximation  for  equation 
(8).  At  higher  velocities  the  higher  velocity  exponent  in  equation  (8) 
results  in  increasingly  higher  loads  than  equation  (7). 

Also  plotted  in  Figure  1  is  the  empirical  relation 

P  =  13.25  V  1,05  (9) 

o 

which  was  derived  in  the  main  body  of  the  report.  This  relation  yields 
much  lower  impact  loads  than  equations  (7)  and  (8).  It  should  be 
pointed  out  that  the  effects  of  work  hardening  and  an  elastic-plastic 
material  response  as  they  pertain  to  static  loading  are  contained  im¬ 
plicitly  in  equation  (9). 


I 


m" 1,1  & . 


Figure  1  Impact  load  calculated  by  three  methods  vs.  impact  velocity  (ZnS, 
3  mm  diameter  glass  spheres). 
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CONCLUSION 

Since  ZnS  exhibits  an  elastic-plastic  response  with  significant  work 
hardening  under  static  loading,  it  was  decided  to  use  equations  (8)  and 
(9)  which  account  for  these  effects,  in  the  analysis  of  the  impact  load 
versus  crack  length  relations.  Equation  (7)  was  not  used  since  it  is 
merely  an  idealized  case  of  equation  (8). 


i 
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Abstract 

Crack  growth  in  alumina  was  simulated  for  the  case  of  small  flaws 
in  larger  surface  crystals  to  determine  the  range  of  crystal  sizes  and 
loading  rates  in  which  single  crystal  fracture  energies  would  be  expected 
to  control  fracture.  This  should  be  expected  only  for  large  crystals  and 


very  high  loading  rates. 


60 


I.  Introduction 

It  has  been  pointed  out  by  Rice  and  co-workers  (1-4)  that,  when 
fractures  originate  at  relatively  small  flaws  in  large  grains  in  poly¬ 
crystalline  ceramics,  fracture  may  be  controlled  by  single  crvstal 
fracture  energies  rather  than  the  fracture  energy  of  the  polycrvstalline 
body.  Comparisons  of  fracture  energies  calculated  from  measured  flaw 
sizes  assuming  negligible  subcritical  crack  growth,  witli  fracture  energies 
measured  on  large  cracks,  were  used  to  support  their  viewpoint.  The  calcu 
lated  fracture  energies  frequently  were  smaller  than  the  measured  values. 
However,  using 


K  =  (2Ey  )1/2  =  Yo  a  1/2  (1) 

Ic  c  r  c 

in  which  is  the  critical  stress  intensity  factor,  E  is  Young's  modulus 

of  elasticity,  is  the  critical  fracture  energy,  Y  is  a  geometrical 
constant  that  accounts  for  the  shape  and  location  of  the  crack,  o  is  the 

r 

fracture  stress,  and  a  is  the  strength  controlling  dimension  of  the  crack 

c 

one  can  see  that,  if  subcritical  crack  growth  occurs,  a^  is  underesti¬ 
mated,  so  that  it  is  uncertain  whether  or  not  the  difference  in  the 
fracture  energies  represents  physical  reality  or  the  result  of  an  incor¬ 
rect  assumption  about  subcritical  crack  growth. 

Recently,  Singh  ££  al_.  (5)  investigated  the  relationship  of  flaw 
size  and  grain  size  to  the  strength  of  ceramics.  They  concluded  that  in 
most  ceramics  an  initial  flaw  contained  within  a  single  large  grain  will 
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initially  propagate  and  then  arrest  after  entering  the  region  of  increased 
fracture  surface  energy.  However,  their  analysis  also  neglected  to 
account  for  subcritical  crack  growth. 

Calculations  of  subcritical  crack  growth  for  glass  and  alumina  in 
water  at  room  temperature  have  shown  that  very  high  loading  rates,  close 
to  those  characteristic  of  Charpy  impact  tests,  are  necessarv  to  limit 
subcritical  crack  growth  to  negligible  values  under  these  conditions  (6). 
Therefore,  the  assumption  of  negligible  subcritical  crack  growth  requires 
further  evaluation.  In  this  investigation  subcritical  crack  growth  calcu¬ 
lations  were  used  to  determine  the  conditions  at  failure  in  specimens  in 
which  the  flaw  at  the  fracture  origin  was  smaller  than  the  particular 
crystal  at  the  origin  and  the  size  of  the  crystal  was  varied  over  a  wide 

—  1  8  —  i 

range.  Loading  rates  were  varied  from  10  to  10  MPas  and  included 
delayed  fracture  tests.  The  procedures  are  described  in  the  next  section. 


62 


I  1  .  Procedures 


Selection  of  data  for  calculations 

In  order  to  provide  a  quantitative  theoretio.il  basis  for  lud-z.  in-.: 
the  importance  of  suber it ical  crack  growth  in  determinat ion  of  fracture 
energies  from  experimental Iv  measured  fracture  stresses  and  flaw  sizes, 
.-.niter  it  ii:al  crack  growth  was  calculated  for  small  flaws  in  larger  grains 
of  various  sizes  under  various  loading  conditions.  i!r>se  calculations 
were  done  by  numerical  integrations  involving,  the  oiu.it  inn 


in  which  V  is  the  crack  velocity,  is  the  stress  intensity  factor,  and 

A  and  n  are  constants  for  a  particular  material  and  set  of  test  conditions 
•k 

Alumina  was  the  only  material  for  which  suitable  data  (values  of  A,  :i , 

and  K,  )  were  known  to  be  available  for  a  single  crvstal  and  a  pel  v- 
1  c 

crvstal  line  body  of  the  same  material,  measured  at  the  same  ti-nnrr.it  uti 
and  humidity.  The  data  of  Wi ederho rn  (7,81  for  sapphire  and  Kv.ins  ( 8 , 9 1 
for  a  coarse  grained  (2S  cm  average  grain  size  (10)1  jm  1  y c rys t a  1  !  ine 
alumina  body  met  these  conditions  in  the  case  of  measurements  at  >-)" 
relative  humidity .  However,  preliminary  calculations  revealed  that 


Despite  the  fact  that  the  choice  of  alumina  was  no  i  -.sar-.  ,  ii  is  subject 
to  certain  disadvantages.  The  principal  d  i  sadvant  a  t  '  it  the  thermal 

expansion  anisotropy  of  corundum  gives  rise  to  1  ■  i  I  :  •«  :  r.  s  idu  ;  1  strrssr 
that  are  not  accounted  for  in  the  analysis. 


f.l 


suhcrit ica 1  crack  growth  was  very  Important  for  this  particular  case 

primarily  because  there  is  such  a  large  difference  between  the  K  at 

which  significant  crack  growth  begins  in  the  sapphire  l.n  MPan'  ) 

!  '  1 

and  the  K  of  the  polvcrvstalline  bodv  (K,  =  p.a  MPam  "l.  Because  the 

ic  •  1c 

above  data  represent  a  particularly  unfavorable  case  for  evaluation  of 

the  hypothesis  involving  the  possibility  of  fracture  of  polvcrvstalline 

bodies  at  single  crystal  fracture  energies,  a  second  set  of  data  was 

calculated  for  a  fine  grained  alumina  body  with  a  subs t ant i a  1  l y  lower  K 

The  body  selected  was  a  96"  alumina.  Values  of  n  and  A  for  Region  I  were 

estimated  using  the  data  of  Freiman,  McKinney,  and  Smith  (2)  that  apparently 

were  taken  at  407  relative  humidity.  The  transition  region  (Region  II) 

which  is  mainly  determined  by  geometry  and  humidity  was  assumed  to  occur 

at  the  same  crack  velocity  as  that  for  Evans'  coarse  grained  alumina. 

K  was  assumed  to  he  the  same  as  that  measured  bv  Bans a  1 ,  Duckworth,  and 
1  c 

Niesz  (11)  for  the  same  alumina  fractured  under  water.  Because  of  the 
high  crack  velocity,  the  presence  of  the  water  should  not  have  a  signif¬ 
icant  ettect  on  the  K  value.  Even  though  the  conditions  used  in  the 

L  e 

three  investigations  from  which  the  data  were  taken  for  the  fine  grained 

alumina  calculations  are  not  entirely  self  consistent,  these  data  have 

enabled  us  to  compare  results  for  two  aluminas  with  differing,  K  values. 

I  c 

* 

The  data  are  given  in  Figure  1  and  Table  I  . 

General  description  of  assumed  c rack  propagation  condi t i ons 

For  tile  growth  of  the  flaw  within  the  crystal  at  the  fracture  origin, 
the  velocities  were  estimated  using  an  approximation  of  the  relationship 

k 

The  symbols  SCI,  SCI1,  SCUT,  and  PCI,  PCI  I,  PCI  I  I  refer  to  Regions  1, 
fl,  and  IJI  for  single  crystals  and  po 1 vervsta 1 s ,  respectively. 


Crack  Velocity  -ms" 


Figure  I  Crack  velocity  vs.  stress  intensity  factor  for 
single  crystal  and  polycrystalline  alumina 
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Table  I 


Values 

of  n,  A 

and  the  end 

points  used  for 

subcritical  c 

rack  growth 

cal  collations 

Crack 

Growth 

Material 

Region 

n 

A 

Find  Point 

Coarse  grained 

alumina  (1) 

SC  I 

112.70 

1.07-  10"9 

Kt  ■  1.02 

SC  II 

0.00 

4.00-  io"c? 

Kt  l.c 

SC  III 

66.  89 

5. 12- 10"19 

Crys  ‘  -il 
Pomviar.v 

PC  I 

33.oo 

A.  98-  10~211 

K-  •:  A  .  12 

PC  II 

0.00 

0.001 

K  il 

rTri 

Firuj  grained 

alumina  (2) 

SC  I 

112.70 

1.07- IQ’9 

KT  •  1 . 02 

SC  II 

o.co 

A.  00-  10" r> 

K,  =  1.6 

SC  [IT 

66. 89 

-IQ 

5.12-10 

C rye i al 
bound a  r.v 

PC  I 

51.  Ao 

_  an 

1.12-  10 

kt  3.3A 

PC  IT 

0. 00 

0.001 

Kt  3.8A 

I  e  ' 

Note  (1)  HYALUMTN/' 

i,  95^  alumina,  25  ji 

m  average  grain 

size 

Note  (2)  ALSIMAG  6lA ,  96 alumina,  5  um  average  (train  size, 
3  M  Company,  cnattanooga,  Tennessee 
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between  crack  velocity  V  and  stress  intensity  depicted  in  the  curve 
(Figure  1)  for  sapphire.  When  the  crack  reached  the  boundary  of  the 
crystal,  the  crack  growth  was  assumed  to  be  governed  by  one  or  the  other 
of  the  curves  for  the  polycrvstalline  materials.  Whether  or  not  the 
assumption  of  an  abrupt  transition  from  single  crystal  to  polycrvstal 
propagation  conditions  is  realistic  may  depend  to  a  great  extent  on  the 
size  of  the  crystal  in  which  the  fracture  originates  at  the  fracture 
origin  and  the  variations  in  the  degree  of  preferred  orientation  near  the 
fracture  origin.  To  take  some  specific  examples,  a  50  urn  grain  at  the 
fracture  origin  in  a  body  with  a  5  urn  average  grain  size  will  probablv  be 
surrounded  by  many  grains  so  that  assumption  of  an  abrupt  transition  may 
be  reasonable.  However,  a  50  pm  grain  at  the  fracture  origin  in  a  body 
with  a  25  pm  average  grain  size  will  be  surrounded  by  a  relatively  smaller 
number  of  grains.  Considering  the  fact  that  fracture  will  originate  at 
the  location  that  is  most  vulnerable  in  terms  of  flaw  size  and  local 
variations  in  one  might  expect  to  observe  some  cases  in  which 

propagation  into  neighboring  grains  is  governed  by  fracture  energies  lower 
than  the  polycrystalline  values. 

Most  of  the  calculations  involved  integration  along  two  or  three 
branches  of  the  single  crystal  crack  velocity  vs.  curve  and  two  branches 
of  the  polycrystal  curve.  Region  II  of  the  polvcrystal  curves  was  assumed 

to  extend  to  K_  where  failure  occurred. 

Ic 

In  addition  to  the  other  assumptions,  it  was  assumed  that  the  ceramics 
were  annealed  before  stressing.  This  assumption  permits  us  to  neglect  the 
contribution  that  large  scale  residual  stresses  resulting  from  machining  or 
other  surface  damage  usually  makes  to  the  stress  intensity  factor  (12,13). 


J 
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Methods  of  numerical  integration 

The  numerical  integrations  were  done  using  a  programmable  printing 
* 

calculator  . 

The  integrations  involve  several  steps.  A  selected  increment  of 
time  (uT)  is  used  to  calculate  the  applied  stress  (n)  using  a  =  o  +  7 A T 
in  which  ao  is  the  stress  applied  in  the  previous  step  and  a  is  the 
stress  rate.  This  is  used  to  calculate  using 

1  /2 

Kj  =  (3) 

in  which  (a)  is  the  crack  size.  Semi-circular  surface  flaws  were  assumed 
and  {  was  taken  equal  to  1.38  which  is  the  value  determined  by  Smith, 
Emery,  and  Kobayashi  (14)  for  the  highest  along  the  boundary  of  a 
semi-circular  flaw.  The  new  crack  velocity  is  calculated  using  Equation 
(2).  The  crack  velocity  is  used  to  increment  the  crack  length  using 
a  -  aQ  VAT  in  which  a^  is  the  crack  length  from  the  previous  step. 

Then,  a  new  stress  is  calculated  and  the  calculation  is  repeated.  The 
computation  stops  when  the  desired  value  of  ,  a,  o,  T,  or  V  is  reached. 

Originally,  the  time  increments  were  selected  in  advance.  In  manv 
cases  the  original  program  was  very  inefficient  because  very  small  time 
increments  were  necessary  to  prevent  overshooting  the  desired  end  point 
or  overflowing  the  calculator.  To  remedy  this  problem  a  loop  in  the 
program  is  used  to  select  an  appropriate  AT  whenever  the  desired  maximum 
increment  of  a  particular  variable  such  as  K  is  exceeded.  This  change 

* 

Hewlett  Packard,  Corvallis,  Oregon,  Model  HP-97. 
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reduces  the  time  necessary  to  integrate  most  curves  with  five  branches 
to  less  than  two  hours. 

Another  precaution  is  that  Equation  (3)  applies  only  for  crack  veloci¬ 
ties  less  than  the  terminal  crack  velocity  in  a  particular  material 

(approximately  0.38  C.  where  C  is  the  longitudinal  wave  velocity  in  the 

o  o 

material).  The  program  compares  the  velocity  calculated  by  Equation  (3) 
with  the  terminal  crack  velocity  and  selects  the  lower  value  for  calcula¬ 
tion  of  the  new  crack  length.  Experience  has  shown  that  this  precaution 
does  not  have  an  appreciable  effect  on  the  final  results  in  most  cases. 

In  these  calculations  the  terminal  crack  velocity  of  sapphire  was  taken 
as  3850  nrs  ^ . 

The  program  was  further  revised  making  it  possible  to  reverse  the 
integration  so  that  the  initial  flaw  size  could  be  estimated  from  the 
calculated  critical  flaw  size.  Comparison  of  calculated  initial  flaw 
sizes,  determined  using  the  fracture  stresses  and  critical  flaw  sizes  of 
fractured  specimens,  with  the  flaw  sizes  measured  on  the  fracture  surfaces, 
provided  added  assurance  that  the  calculations  were  realistic.  These 
calculations  were  made  for  five  delayed  fracture  specimens  and  one  speci¬ 
men  fractured  using  a  constant  loading  rate.  The  specimens  were  96%  A^O^ 
and  were  tested  at  20%  relative  humidity.  Although  there  is  an  error 
involved  in  using  the  values  of  n  and  A,  determined  for  40%  humidity, 
to  calculate  the  flaw  sizes  of  specimens  fractured  at  20%  humidity,  this 
error  is  believed  to  be  small.  The  magnitude  of  the  error  was  estimated 
using  earlier  data  (15)  for  the  strengths  of  96%  alumina  fractured  at 
various  humidities.  The  flaw  sizes  in  specimens  fractured  at  20%  rela¬ 
tive  humidity  are  about  10%  smaller  than  those  cf  specimens  fractured  at 
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40%  relative  humidity  if  comparison  is  made  at  the  same  probabilitv  of 
fracture. 

In  cases  in  which  the  crack  growth  began  in  a  single  crystal  region 
and  the  critical  crack  boundary  is  in  a  polycrystalline  region,  the  size 
of  the  flaw  containing  crystal  must  be  known  and  the  proper  change  in  n 
and  A  made  at  this  boundary.  However,  a  preliminary  calculation  indicates 
that  the  likely  errors  in  estimating  the  crystal  size  change  the  calculated 
initial  flaw  size  only  slightly.  Also,  in  these  cases,  it  is  necessary 
to  limit  the  crack  velocity  in  the  single  crystal  region  to  a  maximum  of 
the  terminal  crack  velocity,  otherwise  the  initial  flaw  sizes  may  have 
negative  values. 


W  .‘jHSpiffB 
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III.  Results 


Constant  stressing  rate 

Assuming  10  pm  flaws  in  surface  crystals  of  various  sizes  and  a 

constant  stressing  rate  of  10  MPas  \  crack  growth  was  calculated  for 

coarse  grained  and  fine  grained  alumina  with  the  results  shown  in  Figures 

2  and  3.  Although  there  is  considerably  less  subcritical  crack  growth 

in  the  finer  grained  material  with  the  lower  K  value,  the  amount  of 

Ic 

subcritical  crack  growth  is  still  substantial  with  the  critical  crack 
lengths  ranging  from  400-1400%  of  the  original  flaw  size  depending  on 
the  size  of  the  crystal  containing  the  original  flaw.  In  the  coarser 
grained  material  the  critical  crack  lengths  ranged  from  600-2800%  of  the 
original  flaw  size. 

At  a  stressing  rate  of  10  MPas  \  crystal  sizes  greater  than  75  urn 
in  the  finer  grained  material  and  100  pm  in  the  coarser  grained  material 
were  necessary  before  the  crack  would  "go  critical"  with  respect  to  the 
polycrystalline  K  in  the  flaw  containing  crystal,  so  that  the  crack 
was  not  slowed  down  noticeably  at  the  crystal  boundary. 

Variable  stressing  rate 

The  results  of  the  constant  stressing  rate  calculations  were  used 
to  select  crystal  sizes  for  variable  stressing  rate  calculations  for  simu¬ 
lation  of  crack  propagation  at  various  stressing  rates.  The  chosen  crystal 
sizes  were  50  pm  and  100  pm.  The  stressing  rates  were  chosen  to  indicate 


Crack  Length-pm 


Loading  Time -seconds 


figure  2  Crack  length  vs.  time  for  10pm  flaws  in  crystals  of 

various  sizes  in  a  coarse  grained  alumina  (10  MPas"1 
stressing  rate,  50%  relative  humidity). 


Crack  Length -pm 


Loading  Time -seconds 


Figure  3  Crack  length  vs.  time  for  10pm  flaws  in  crystals  of 
various  sizes  in  a  fine  grained  alumina  (lOMPas-1 
stressing  rate,  50%  relative  humidity). 
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the  range  of  stressing  rates  for  which  the  crack  behaved  as  though  it  was 
growing  mainly  in  a  medium  with  polycrystalline  fracture  properties  and 
the  range  of  stressing  rates  in  which  it  behaved  as  though  it  was  growing 
mainly  in  a  medium  with  single  crystal  fracture  properties.  The  results 
are  given  in  Figures  4-7. 

In  the  fine  grained  alumina  with  a  50  pm  crystal  at  the  fracture 

4-1 

origin,  stressing  rates  greater  than  10  MPas  are  necessary  to  suppress 
subcritical  crack  growth  sufficiently  so  that  reached  K  before  the 

crack  reached  the  boundary  of  the  50  pm  crystal.  With  a  100  pm  crystal  at 

? 

the  fracture  origin,  stressing  rates  greater  than  10“  MPas  were  required 

In  this  case,  even  at  lower  stressing  rates,  the  rate  of  crack  growth  at 

the  100  pm  crystal  boundary  is  great  enough  so  that  the  increases  in 

are  caused  mainly  by  increases  in  crack  length  so  that  there  is  little  or 

no  inflection  of  the  a  vs.  curves  as  they  cross  the  boundary  from 

single  crystal  to  polycrystalline  fracture  properties. 

In  the  coarse  grained  alumina  with  a  50  pm  crystal  at  the  fracture 

origin,  the  higher  K  of  the  polycrvstalline  material  results  in  a 

greater  range  of  stress  intensity  factors  in  which  the  crack  is  essentiall 

arrested  at  the  crystal  boundary.  Even  at  a  stressing  rate  of  10^  MPas  ^ 

this  range  is  substantial.  Therefore,  it  appears  that  no  reasonable 

stressing  rate  is  sufficient  to  cause  KT  to  reach  KT  before  the  crack 

I  Ic 

reaches  the  boundary  of  the  crystal.  With  a  100  pm  crystal  at  the  frac- 

4  -1 

ture  origin,  stressing  rates  greater  than  10  MPas  were  required. 


Crack  Length-pm 


Figure  4  Crack  length  vs.  stress  intensity  factor  for  various 
stressing  rates  in  fine  grained  alumina;  10 pm 
flaw  in  a  50pm  crystal. 
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Figure  5  Crack  length  vs.  stress  intensity  factor  for  various 
stressing  rates  in  fine  grained  alumina;  10  pm 
flaw  in  a  100pm  crystal. 
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The  labels  on  the  curves  are 
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Stress  Intensity  Factor -MPa m,/2 


Figure  6  Crack  length  vs.  stress  intensity  factor  for  various 
stressing  rates  in  coarse  grained  alumina;  10pm  flaw 
in  a  50 pm  crystal. 


Stress  Intensity  Factor  -  MPa  ml/2 


Figure  7  Crack  length  vs.  stress  intensity  factor  for  various 
stressing  rates  in  coarse  grained  alumina-,  10pm 
flaw  in  a  100pm  crystal. 
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Tin*  results  are  given  in  Figure  H  and  are  compared  wi  th  results  .  ■ :  sir.;',  ar 
calculat  ions  in  which  the  flaw  was  assured  t  >  extend,  te  the  heundurv  of  the 
ervsta]  when  the  loading  began .  The  results  shew  that  ever  a  wide  range 
of  times  to  failure  there  is  no  significant  difference  between  the  tires 
to  failure  for  the  10  ;.m  flaws  in  the  50  jis  ervstals  and  the  50  ..p-  flaws 
extending  to  the  ervsta 1  boundaries.  evidently,  the  tine  required  f.*r 
the  flaw  to  grow  to  the  crystal  boundary  is  negligible  compared  with  tin- 
time  to  failure. 

The  delayed  fracture  curve  for  the  fine  grained  alumina  has  a  lower 
slope  than  that  of  the  coarse  grained  alumina.  As  a  result,  at  stresses 
less  than  about  225  MPa,  the  times  to  failure  of  the  fine  grained  alumina 
are  greater  than  those  of  the  coarse  grained  bodv  despite  the  much  higher 
K  of  the  coarse  grained  body.  The  shaded  area  in  Figure  8  represents 
results  of  delayed  fracture  experiments  using  the  fine  grained  96?  alumina. 
The  upper  limit  is  for  the  laboratory  atmosphere  and  the  lower  limit  is 
for  a  water  environment.  The  fact  that  the  calculated  curve  passes  through 
this  shaded  area  with  an  intermediate  slope  shows  that  the  calculations 
are  realistic.  The  96Z  alumina  used  for  these  experiments  was  in  the 
as-fired  condition.  This  condition  is  consistent  with  the  assumption  of 
negligible  large  scale  residual  stresses  at  the  surface  flaws. 


Calculated  results  for : 


Log  Time  to  Failure -seconds 

Figure  8  Fracture  stress  vs.  log  time  for  fine  and  coarse  grained  aluminas 

assuming  either  a  10 pm  flaw  in  a  50pm  crystal  or  a  50pm  flaw  in  the 
polycrystalline  body. 
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Comparison  of  calculated  and  measured  values  of  initial  flaw  size 

Critical  flaw  sizes  in  96 7.  alumina  specimens,  fractured  in  an  earlier 
program  (16,17)  were  calculated  using  Equation  (1)  in  which  semi-circular 
surface  flaws  were  assumed  and  Y  was  taken  as  1.27.  I'sing  these  critical 
flaw  sizes,  the  integrations  were  reversed  to  yield  estimates  of  the 
initial  flaw  sizes.  These  flaw  sizes  were  compared  with  the  sizes  of  the 
flaws  at  the  fracture  origins,  observed  in  scanning  electron  micrographs 
of  the  fracture  surfaces.  The  results  obtained  for  several  delaved  frac¬ 
ture  specimens  are  given  in  Table  II.  Reasonable  agreement  between  the 
calculated  and  measured  values  of  initial  flaw  size  was  observed.  The 
flaws  in  specimens  D-3,  D-ll,  and  D-12  were  stepped  flaws  whereas  D-6 
and  D-7  were  large  crystals  of  30  um  and  25  pm,  respectively.  In  D-7 
no  flaw  was  visible  in  the  large  crystal  so  the  entire  crystal  was 
assumed  to  be  the  flaw.  In  D-6,  a  small  flaw  was  clearly  visible  at 
the  surface  of  the  crystal.  In  this  case  the  calculation  indicated  that, 

at  the  crystal  boundary,  the  crack  growth  rate  drops  from  3850  m-s  ^  to 
_8 

7. 5*10  m-s  .  The  results  obtained  for  specimen  D-6  represent  the  best 
test  of  these  crack  growth  estimates.  In  this  case  the  original  flaw  size 
was  5  um  in  a  30  um  crystal.  The  calculated  critical  flaw  size  was  51.1  n m, 
1020%  of  the  original  flaw  size,  which  is  roughly  in  the  middle  of  the 
range  expected  based  on  the  constant  stressing  rate  calculations. 

In  the  case  of  a  specimen  fractured  using  a  constant  stressing  rate, 
the  calculated  initial  flaw  size  was  28.6  um.  This  flaw  size  is  in  reason¬ 
able  agreement  with  the  size  of  the  large  crystal  (20  um)  observed  at  the 
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fracture  origin.  The  stressing  rate  used  in  this  calculation  was  based 
on  the  assumption  that  failure  occurred  in  one  minute. 

The  results  in  this  section  show  that  the  calculations  yield  reason¬ 
able  estimates  of  the  initial  flaw  sizes  and  that  it  is  unlikely  that 
these  fractures  are  controlled  by  single  crystal  fracture  energies  even 
in  the  case  in  which  fracture  originated  at  a  5  pm  flaw  in  a  30  pm  crystal. 
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IV.  Discussion  and  Conclusions 

A  basic  assumption  on  which  these  calculations  are  based  is  that 
when  the  crack  passes  the  crystal  boundary,  the  crack  velocity  will 
immediately  or  almost  immediately  change  to  that  expected  based  on  the 
polycrystalline  V  vs.  K  curve.  The  energy  acquired  by  the  crack  while 
it  is  in  the  crystal,  in  excess  of  that  required  for  crack  propagation 
in  the  crystal,  depends  on  the  elastic  energy  input  to  the  crack  and  is 
partitioned  among  the  stress  wave  energy  (5),  the  excess  surface  energy 
absorbed  by  K  dependent  mechanisms,  and  the  kinetic  energy.  Clark  and 
Irwin  (18)  refer  to  running  cracks  as  low-inertia  disturbances  that 
cannot  be  driven  faster  than  about  half  of  the  elastic  shear  wave  velocity. 
Therefore,  it  is  reasonable  to  assume  that  the  kinetic  energy  is  small  and 
has  a  definite  upper  limit  depending  on  the  shear  wave  velocity.  Conser¬ 
vation  of  energy  requires  that,  after  the  crack  crosses  the  crystal 
boundary  into  the  polycrystalline  region,  the  energy  required  to  make  up 
the  difference  between  the  elastic  energy  input  and  the  fracture  energy 
absorbed  in  this  region  which  requires  much  greater  fracture  energy  for 
propagation,  must  be  drawn  from  the  kinetic  energy. 

Quantitative  methods  are  not  available  to  permit  determination  of  the 
kinetic  energy.  However,  it  is  possible  to  estimate  an  upper  limit  to 
crack  propagation  before  "arrest"  in  the  polycrystalline  region  by 
equating  the  excess  elastic  energy  supplied  to  the  crack  tip  while  it  is 
in  the  crystal  to  the  excess  fracture  energy  absorbed  in  the  polycrystalline 
region.  This  calculation  requires  integration  over  the  area  of  crack  growth 
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in  each  case.  This  was  done  for  the  case  of  a  crack  that  went  critical 
at  a  crack  radius  of  25  ijid  in  a  50  pm  crvstal.  At  a  radius  of  56  urn, 
the  excess  elastic  energv  had  been  entirely  absorbed  as  fracture  energy. 
These  results  show  that  it  is  reasonable  to  assume  that  the  crack  slows 
down  to  the  crack  velocity  expected  based  on  the  polycrvstalline  V  vs. 
curve  as  the  crack  passes  the  crystal  boundary. 

The  results  of  the  crack  growth  calculations  show  that,  in  alumina 

ceramics  with  small  flaws  in  large  grains,  there  is  a  range  of  flaw  sizes, 

grain  sizes,  and  stressing  rates  in  which  the  crack  growth  becomes  critical 

with  respect  to  the  polvcrystalline  KT  before  the  crack  reaches  the 

lc 

boundary  of  the  grain  in  which  the  fracture  originates.  However,  in  the 
cases  examined  thus  far,  substantial  subcritical  crack  growth  occurred, 
significantly  changing  the  critical  flaw  size  and  fracture  stress.  This 
was  true,  even  at  very  high  stressing  rates.  This  finding  is  in  basic 
disagreement  vith  the  assumption  of  Rice  and  co-workers  that  subcritical 
crack  growth  can  be  neglected  so  that  fracture  energies  can  be  calculated 
from  initial  flaw  sizes  at  fracture  origins  and  the  fracture  stresses. 

The  effects  of  neglecting  subcritical  crack  growth  on  fracture  energies 
are  indicated  in  Figure  9.  The  calculated  fracture  energies  are  called 
apparent  fracture  energies.  As  indicated  by  the  solid  line  in  the  figure, 
the  apparent  fracture  energies,  calculated  based  on  the  original  flaw 
size  and  assuming  that  the  crack  propagation  is  controlled  by  the  poly¬ 
crystalline  properties  (A  and  n) ,  are  lower  than  polycrystalline  fracture 
energies  but  not  drastically  so.  However,  in  cases  in  which  subcritical 
crack  growth  occurs  in  relatively  large  single  crystals  at  the  fracture 
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origins  (indicated  by  the  dashed  lines  for  50  and  100  ,.m  crvstals)  there 
is  a  drastic  reduction  in  the  apparent  fracture  energy.  This  reduction 
is  caused  by  the  extensive  subcritical  crack  growth  that  occurs  in  these 
cases  which  reduces  the  calculated  fracture  stresses.  These  low  fracture 
stresses,  when  combined  with  the  underestimated  flaw  sizes  that  result 
from  using  the  original  flaw  sizes  instead  of  the  critical  flaw  sizes, 
yield  substantial  underestimates  of  the  fracture  energies. 

The  present  method  of  calculation  can  be  used  to  investigate  the 
probability  that  fractures  in  alumina  will  be  controlled  by  the  single 
crystal  in  practical  applications.  It  seems  clear  that  most  practical 

applications  involve  static  loads  or  low  stressing  rates.  As  indicated  bv 
the  delayed  fracture  simulations,  it  is  unlikely  that,  under  these  con¬ 
ditions,  the  fractures  will  be  controlled  bv  the  single  crystal  .  In 

fractures  occurring  at  high  loading  rates  (ceramic  armor,  impacts  of 
pellets  on  high  intensity  lighting  envelopes  or  electrical  insulation, 
interrupted  cutting  with  ceramic  cutting  tools,  etc.),  control  by  the 
single  crystal  is  more  likely.  Kirchner,  Gruver,  and  Softer  (16,17) 

characterized  flaws  at  fracture  origins  in  967.  alumina  and  H.  P.  alumina 
for  various  loading  rates  and  temperatures.  Depending  on  the  conditions, 
the  percentage  of  fractures  originating  at  large  crystals  ranged  from  0 
to  657.  but  frequently  was  about  20-307..  These  crystals  were  less  than 
50  jim  in  size  with  possible  rare  exceptions.  Therefore,  present  evidence 
indicates  that  in  fractures  originating  at  large  crystals,  the  poly¬ 
crystalline  Kj  will  control  fracture  in  most  cases  and  that  the  single 
crystal  will  control  fracture  only  in  the  case  of  the  largest 
crystals  and  the  highest  loading  rates. 
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Further  thought  should  be  given  to  several  other  effects  including: 

1.  The  effect  of  residual  stresses  at  the  machining  flaws. 

2.  The  possible  effect  of  localized  variations  in  preferred  orien¬ 
tation  on  localized  KT  . 

Ic 

Residual  stresses  can  have  a  profound  effect  on  the  K^.  values  actually 

present  at  crack  tips.  Petrovic  £t  al_.  (12)  found  that  residual  stresses 

at  Vickers  indentations  in  H.P.  Si.N,  reduced  measured  KT  to  about  3.5 

2  4  Ic 

1/2  1/2 

MPam  from  the  value  of  4.7  MPam  otherwise  expected.  Kirchner  and 

Isaacson  (13)  found  a  similar  reduction  in  the  case  of  single  point  diamond 

grinding  damage  where  measured  values  were  reduced  to  about  2.8 

1/2 

MPam  .  In  addition,  stresses  caused  by  thermal  expansion  anisotropy 
may  have  a  more  localized  effect.  If  these  effects  are  as  large  in 
alumina  as  in  H.P.  Si^N^,  it  seems  possible  that  in  some  cases  the  original 
flaws  might  propagate  to  the  grain  boundary  under  the  influence  of  the 
residual  stresses  alone. 

The  presence  of  preferred  orientation  will  probably  result  in  a 
gradual  transition  from  single  crystal  to  polycrystalline  fracture  energies 
near  the  crystal  boundary.  In  evaluating  this  possibility  it  is  well  to 
remember  that,  in  effect,  the  fracture  seeks  out  the  most  vulnerable  frac¬ 
ture  origin  in  terms  of  flaw  size  and  severity,  susceptibility  to 
subcritical  crack  growth  and  local  .  Thus,  if  there  are  two  similar 
flaws  in  two  similar  crystals  but  one  is  surrounded  by  a  region  of  pre¬ 
ferred  orientation  that  is  favorable  for  subcritical  crack  growth  and  the 
other  is  surrounded  by  randomly  oriented  crystals  the  resulting  fracture 
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